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Abstract 



A Forward Plug Calorimeter (FPC) for the ZEUS detector at HERA 
has been built as a shashlik lead-scintillator calorimeter with wave 
length shifter fiber readout. Before installation it was tested and cali- 
brated using the X5 test beam facility of the SPS accelerator at CERN. 
Electron, muon and pion beams in the momentum range of 10 to 100 
GeV/c were used. Results of these measurements are presented as well 
as a calibration monitoring system based on a 60 Co source. 
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1 Introduction 



The ZEUS |jj collaboration has installed a Forward Plug Calorimeter (FPC) 
around the beam line (see Fig. ||) to extend the calorimetric coverage in 
pseudorapidity [] from r/ < 4.0 to 77 < 5.0. This vastly increases the physics 
potential for diffraction in deep inelastic scattering. The mass range over 
which the dissociated photon system can be studied is enhanced by a factor 
of about two ||]. 

This document describes the construction and the results of the beam test 
of the FPC performed at CERN in September 1997 with electrons, muons 
and pions covering an energy range from 10 to 100 GeV. Prior to this, 
a prototype corresponding to one half of the electromagnetic section was 
tested at DESY with electron beams from 1 to 6 GeV Data from the 
beam test at DESY will also be presented where appropriate. 

2 Description of the detector 

The FPC is a lead-scintillator sandwich calorimeter read out by wave length 
shifter (WLS) fibers and photomultipliers (PMT). This concept has been 
investigated in j|, @j. It has been installed in the 20 x 20 cm 2 beam 
hole of the forward uranium-scintillator calorimeter (FCAL) of the ZEUS 
detector at HERA (see Fig. |). A front view of the FPC is shown in 
Fig. g. The FPC has a 63 mm diameter central hole to accommodate the 
HERA beampipe. 

The active part of the FPC has outer dimensions of 192 x 192 x 1080 mm 3 . 
The FPC is mechanically subdivided into two identical half modules. They 
are attached to the bottom and top half parts of the innermost FCAL mod- 
ule Q. The FCAL halves and therefore also the FPC halves are moved apart 
for beam injection in order to reduce the radiation dose on the calorimeters. 

In the FPC, lead plates of 15 mm thickness alternate with scintillator layers 
of 2.6 mm. The WLS fibers have 1.2 mm diameter and pass through 1.4 mm 
diameter holes in the lead and scintillator layers. The holes are located on 
a 12 mm step square grid. The FPC has 232 holes of this type. In addition, 

1 Pseudorapidity is defined as r\ — — In (tan(#/2)), where 6 stands for the angle between 
the particle trajectory and the forward proton beam direction. 
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Figure 1: Side view of FPC integrated in the ZEUS Calorimeter. 

there are 4 brass tubes (1.4 mm inner diameter) which are used to guide a 
60 Co source, placed on the tip of a long steel wire, for the monitoring of the 
calibration of the FPC cells. Taking into account the WLS fibers and the 
fiber holes in the lead plates, the effective lead to plastic ratio by volume 
is 5.2:1. Based on the results from a lead-scintillator calorimeter of similar 
composition the FPC is expected to provide equal response to electrons 
and hadrons (compensating calorimeter, e/h = 1). 

With the layer structure chosen, the FPC has approximately the same 
radiation length Xq and nuclear absorption length A as the FCAL, viz. 
Xo(FPC) = 0.68 cm and A (FCAL) = 0.74 cm; A(FPC) = 20 cm and 
A(FCAL) = 21.0 cm. This minimizes the fluctuations in the energy mea- 
surement. 

The FPC is subdivided longitudinally into an electromagnetic (EMC) and 
a hadronic (HAC) section which are read out separately (see Figs. || and [|). 
The electromagnetic section consists of 10 layers of lead and scintillator 
corresponding to 26.5 Xq and 0.9 A. The hadronic section of the FPC consists 
of 50 layers and represents 4.5 A leading to a total for the FPC of 5.4 A (see 
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Table |). 



The scintillator layers consist of tiles and form cells which are read out 
individually (Fig. ||). The cell cross sections are 24 x 24 mm 2 in the EMC, 
commensurate with the Moliere radius, and 48 x 48 mm 2 in the HAC section. 
The 8 (4) innermost cells in EMC (HAC) surrounding the beam hole follow 
the circular shape given by the beam hole. 

The polystyrene based scintillator SCSN81T2 from Kuraray was used, since 
it was found to be the best choice in terms of light yield and radiation stabil- 
ity The scintillator tiles were wrapped with 0.2 mm thick tyvek paper 
in order to improve the light collection efficiency and avoid light coupling 
between neighboring cells. 

The WLS fibers of the EMC are connected to clear fibers (1.4 m length) 
which transport the light to the PMTs placed behind the FPC (see Fig. [|). 
The clear fibers are bent by 180° at the front of the EMC, and are guided to 
the PMTs in the rear by two ducts situated on either side of the FPC. The 
WLS fibers of the HAC section transport the light directly to the PMTs. All 
4 (16) fibers corresponding to an EMC (HAC) cell are connected to the same 
PMT through a light-mixer bar. An additional fiber is connected to each 
PMT through the same light-mixer bar in order to inject LED and laser light 
pulses for monitoring the stability of the PMTs and readout electronics. On 
the other end of the WLS fibers a reflective aluminized mylar foil is placed 
to avoid light losses (see Fig. |3|). The total number of readout channels is 
(EMC + HAC): 60 + 16 = 76 (see Table ||). 

For the WLS fibers the material Y11200(dc) from Kuraray has been chosen. 
It has a polystyrene core and a double cladding (PMMA and fluorinated 
PMMA) which produces a substantial increase in light output compared to 
a single clad fiber. Its absorption spectrum matches best the SCSN81T2 
emission spectrum. 

The PMT chosen for the FPC is the Hamamatsu R5600U, a tube with metal 
channel dynodes. It is relatively insensitive to magnetic fields and its small 
dimensions (16 x 16 x 16 mm 3 ) are well suited for the limited space available 
around the hole of the FCAL. Tests done prior to the beam tests showed 
that it fulfills well the constraints on dynamic ranges placed by the required 
FPC performance in the HERA environment ||]. 

The PMTs are placed in four separate iron blocks whose shape and trans- 
verse size is similar to that of the lead plates. The four structures correspond 
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Figure 2: Front view of the FPC. The readout cells and the position of WLS 
fibers are shown. Notice that there is one hadronic readout cell behind 2x2 
electromagnetic cells, except for the cells near the beam hole. 
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Figure 3: Schematic side view of the FPC. The horizontal lines indicate the 
positions of the WLS fibers and of the tubes for the 60 Co source. 
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0.9 + 4.5 = 5.4 A 



Table 1: Summary of FPC parameters. 



to the EMC and HAC sections of both halves. These blocks are placed be- 
hind the hadronic section, first the one containing the HAC PMTs and then 
that for the EMC (see Fig. |). The EMC PMTs are arranged in their sup- 
port frames in such a way that they do not line up with their respective 
calorimeter cell. 

After the installation in ZEUS, the PMTs are connected to a high voltage 
system based on a Cockroft- Walton generator M. 

3 CERN setup 

The FPC has been tested in the X5 beam of the CERN SPS. It has been 
installed between modules of the FCAL prototype, see Fig. [|. These are 
uranium-scintillator calorimeter units similar to the ones implemented at 
ZEUS (see |l(|]). The uranium plates have a thickness of 3.3 mm and the 
scintillator layers are 2.6 mm thick. The FCAL units are divided longitu- 
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Figure 4: Front view of the FPC at the test beam site, surrounded by three 
FCAL modules. 



dinally into 3 sections, one electromagnetic, 26 Xq, and 2 hadronic sections, 
3 A each. The electromagnetic (hadronic) sections are divided transversally 
into 5 x 20 cm 2 cells (20 x 20 cm 2 cells). 

The combined FPC + FCAL prototype setup has been placed on a structure, 
which could be moved in the x and y directions, in order to vary the impact 
point of the particle beam at the calorimeter. The beam line defines the 
z-direction. 
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Figure 5: Top view sketch of the beam test setup. 
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A Delay Line Wire Chamber (DLWC) in front of the FPC has given precise 
information on the charged particle position (see Fig. ||). The relative po- 
sition of the FPC detector with respect to the DLWC coordinates has been 
established to within about one millimeter. Cuts on DLWC coordinates 
were used in the offline analysis in order to select particles hitting a given 
region of the detector. 

A defocused beam was used to obtain events covering a wide region of the 
FPC. The trigger was defined by requiring signals in both scintillator coun- 
ters (61 and 62) placed in front of the FPC. Behind the FPC an additional 
counter (65) was installed mainly for offline identification of muons. 



4 Simulation of the FPC measurements performed 
at CERN 



For the analysis of the CERN-test data the FPC has been simulated by 
Monte Carlo (MC) using the GEANT 3.21 package 0. The FPC has 
been implemented as a sampling calorimeter consisting of scintillator and 
absorber layers. The density of the absorber layers has been reduced in order 
to account for the tyvek layer and for the holes left for the WLS fibers. The 
fibers were not included in the simulation. The light attenuation in the WLS 
fibers has been taken into account according to the measurement described 
in section 5.4.3] . Photostatistics, noise and cross talk have also been included 



according to measurements results (see section 5.3). In the simulation the 
FPC has been positioned between FCAL modules as in the CERN-test setup 
described in the previous section (see also Fig. Q). 

The accuracy of the simulated response to electrons and muons has been 
estimated to be about 3%. For pions the accuracy of the simulated energy 
signals, which have been determined with the hadronic package GHEISHA 
[13, is about ±10%. 



5 Calibration 



For the calibration of the FPC the following scheme has been chosen: 



• calibration of EMC cells using electron test beam data, 
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• calibration of HAC section using muon test beam data, 



• monitoring of the calibration constants using the response of each FPC 
channel to irradiation with a 60 Co source. 



5.1 Calibration of the EMC section 

A data set of 60 GeV test beam electrons distributed over the complete FPC 
surface has been used for the calibration of the EMC cells. The calibration 
constants have been adjusted such that the mean of the summed signal of all 
EMC cells is independent of the point of incidence. The absolute calibration 
has been adjusted to the beam energy. For the edge cells the beam energy 
is reduced by a correction factor, estimated by MC, taking into account the 
energy leaking out of the FPC. 



5.2 Calibration of the HAC section 



Since the energy of pions is not fully contained in the FPC and the amount 
of leakage is not reliably described by the MC, muons have been used for the 
calibration of the HAC cells. A Landau function convoluted with a gaussian 
function, to take into account photostatistics and electronic noise, has been 
fitted to the muon signals. The calibration constants have been computed in 
order to adjust the peak value obtained from the fit, to the value predicted 
by MC. Since the distribution of the muon signals does not have a gaussian 
shape, it is crucial that it is well described by MC. A good overall agreement 
is observed for the HAC section after calibrating with muons, see Fig. ||. 

In order to measure hadronic energy, the following two corrections have been 
applied: 



The ratio between the visible and the total energy deposited in the 
FPC (sampling fraction) is different for an incident pion and an in- 
cident muon. Assuming that the FPC is a compensating calorimeter 
(e/h = 1), the sampling fraction obtained with electrons (see sec- 
tion |5.1|) can also be used for pions. 



The effect of the light attenuation in the WLS fibers (see section 5.4.3| ) 



is different for muons and pions due to their different longitudinal 
energy deposition. Whereas muons deposit energy uniformly along the 



S 




Figure 6: The reconstructed energy measured in the EMC and HAC sections 
of the FPC for 100 GeV muons is shown for data and MC. 



z-direction, pions deposit their energy predominantly at the beginning 
of the HAC region where the effect of the light attenuation is strongest. 
In order to correct for the effect, the energy scale of the HAC section 
has been raised by a factor of 1.33. 

The accuracy of the absolute calibration thus obtained is estimated to be 
5%. The error is dominated by the uncertainty of the MC prediction. 



5.3 Cross talk 

A cross talk effect has been observed between physically adjacent cells of the 
FPC in the analysis of the CERN-test data. No cross talk between PMTs 
or front-end electronic channels has been found. The observed cross talk 
has therefore been attributed to light being produced in one cell, crossing 
the tyvek barrier between scintillator tiles and entering the WLS fibers of a 
neighboring cell. It is known that tyvek is not totally light tight. 

The cross talk has been measured using the muon sample. For all pairs of 
neighboring cells the ratio between the signal of the cell that has been hit 
and the signal in the adjacent cell has been determined. The mean value of 
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EMC 


HAC 


central cell 


73% 


81% 


each of 4 direct neighbors 


5.8% 


4.3% 


each of 4 diagonal neighbors 


1% 


0.4% 



Table 2: Fraction of light collected in the central cell and in the adjacent 
ones. 



this distribution has been taken as the amount of cross talk. 

For light produced in the scintillator of a given cell (denoted as the central 
cell) the fraction of light that has been collected in this and the adjacent 
cells is shown in Table [2[ 

Note that the method applied to calibrate the EMC (HAC), which has been 
described in previous sections, is not sensitive to cross talk, since it uses the 
sum of the signals over all cells. 

In order to describe the lateral shower width, the measured cross talk has 
been implemented in the MC. 

5.4 60 Co monitoring system 

A monitor system using a 60 Co source, similar to the one used for the ZEUS 
uranium calorimeter fl3|| , has been developed for the FPC. It allows the 
detection of changes in the performance of the scintillator tiles and the WLS 
fibers as well as drifts in the gain of the PMTs. By measuring the ratio of 
response to 60 Co and beam particles the absolute and cell-to-cell calibration 
constants can be transported from the test beam to ZEUS and the stability 
of the calibration can be monitored. 

5.4.1 60 Co setup 

A 1 mCi pointlike 60 Co source is attached to the tip of a steel wire. For 
safety reasons the wire and the source are enclosed inside a stainless-steel 
tube. This source wire can be inserted into brass tubes of 1.4 mm inner 
diameter inside the FPC. The brass tubes run parallel to the WLS fibers 
at the center of each FPC quarter (see Fig. ||). The 60 Co source irradiates 



10 



Co - source (1 .2 MeV y -source) 



e:\ic 



_LLL 



TTTT 
I I I I 
T 



i 



WLS -fibers 



W 

scintillators 



Figure 7: Schematic drawing of the FPC with inserted 60 Co source 

the scintillators of the FPC with 1.173 MeV and 1.332 MeV photons (see 
Fig. §. 

A PC-controlled motor moves the 60 Co source wire in steps of 0.6 mm 
through the FPC. At each step the PMT currents are read out across 2 Mfi 
resistors by an integrating, voltage sensitive ADC. The measurement is re- 
peated 500 times and the mean and root mean square (RMS) values are 
recorded. During a complete scan the source fully traverses the FPC (from 
EMC to HAC) and the 15 EMC and 4 HAC channels of the corresponding 
FPC quarter are read out. 

The signal as a function of position is shown in Fig. |8|. The 10 (50) peaks in 
the EMC (HAC) part correspond to positions where the source is next to a 
scintillator. When the source is inside a lead layer the photons are partially 
shielded and the signal drops. The point where the source exits the EMC 
section and enters the HAC section is clearly seen at z ~ 190 mm. 

The maxima of the signals from the individual scintillator tiles fluctuate. 
This is assumed to be due to differences in the scintillator- WLS fiber light 
coupling and to differences in the tyvek wrapping and scintillator machining. 
The overall increase of the signal in the HAC section with increasing z is 
due to light attenuation in the WLS fibers: the fibers are read out in the 



direction of positive z (see section 5.4.3) 
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Figure 8: The measured 60 Co signal in the EMC and HAC sections as 
function of the source position. The vertical lines indicate the region 
integration used to obtain the value of COMEAN. 
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Figure 9: T/ie COMEAN values from different m Co scans, listed by dates, 
are shown for HAC cell 62 and for three EMC cells at different distances to 
the cobalt source. According to the classification in Fig. cell 7 belongs to 
group A, cell 3 to B and cell 4 to C. 

5.4.2 Monitoring of the stability of the signal 

In order to study the stability of the FPC response, the signal induced by the 
60 Co source is integrated between the bounds shown in Fig. ||. The result of 
the integration is called COMEAN. Figure | shows COMEAN for 8 different 
60 Co scans spanning the time from the CERN-test period in September 1997 
until the installation of the FPC in ZEUS and the start of luminosity runs in 
August 1998. Up to 10% deviations can be observed after the transportation 
of the FPC to DESY (05.11.97) and after installation of the FPC into the 
ZEUS environment (25.04.98). 
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Figure 10: Front view of a quarter of the FPC. The EMC cells are classified 
in three groups depending on their distance to the 60 Co source. 

The different signal heights observed for the different cells in Fig. ^ is due 
to their distances from the cobalt source. The EMC cells can be subdivided 
into 3 groups (A,B,C) depending on their distance from the 60 Co source (see 
Fig. |i~0| ). The distance between the 60 Co source and HAC cells is the same 
for all HAC cells. 



5.4.3 Measurement of the light attenuation in WLS fibers 

The change of the peak height with the source position z in the HAC section 
(see Fig. |8|) allows the measurement of the light attenuation length in the 
WLS fibers directly. Since the peak height is also influenced by the scintil- 
lator quality and the light coupling into the WLS fiber, average peak values 
of all channels have been used (see Fig. [IT| ). 

The data are fitted by using the function: 

f(z) = Pi • (e-C*)-*)/^ + P 3 ■ e -(z-zo+2-L)/P^ (-Q 

where L = 89.1 cm is the length of the HAC (distance from the reflective 
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Figure 11: The peak height averaged over all cells is plotted as a function of 
peak position z. The attenuation length in the WLS fibers is determined with 
a two exponential fit. The second exponential accounts for the light reflected 
at the end of the fiber. 

mylar foil to the last scintillator layer on the back) and zq = 108.0 cm is the 
z-position where the fiber exits the FPC. In Eq. [I] the second exponential 
accounts for the light reflected at the mylar foil which separates the EMC and 
HAC sections; Pi, Pi and P3 are free parameters in the fit. The attenuation 
length measured using this method is P2 = 107 cm. The weight of the 
second exponential function is P3 = 0.80, which indicates that 80% of the 
light is reflected by the mylar foil. 

6 Effect of WLS fibers on the signal 
6.1 Signal uniformity 

The uniformity of the response across the surface of a FPC cell is affected by 
the presence of the WLS fibers. Light produced in scintillator regions close 
to a WLS fiber is collected with a better efficiency than in other regions as 
discussed in [||]. 

The effect of the enhanced FPC response to particles incident close to a 
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Figure 12: The energy measured by the FPC, Eppc, for 100 GeV electrons 
is plotted as a function of the position of incidence. Two 3 mm wide bands 
(measured in y-direction) of position of incidence have been selected, one 
crossing the WLS fibers (solid dots), the other avoiding the fibers (open 
squares). In the upper plot the FPC surface has been perpendicular to the 
beam axis while in the lower plot it has been tilted by 1.7°. 
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WLS fiber has been studied using high energy electrons for which the energy 
resolution [ue/E) is best. In Fig. 12 the FPC response to 100 GeV electrons 
is shown as a function of the point of incidence. The average response is 
enhanced by ~ 15% in regions close to the WLS fibers. This effect is reduced 
for runs where the FPC was tilted by 1.7° (compare Fig. ^ top and bottom) 
corresponding to the mean angle of incidence of beam particles coming from 
the ZEUS interaction point. 



6.2 Effect on the energy resolution 

The enhancement of the signal in the fiber regions also affects the energy 
resolution of the FPC as seen in Fig. |lj|. The distribution of the total 
reconstructed energy, Eppc-, shows a tail to higher energies. Figure |^ right 
shows the positions of electrons selected from this tail, as measured by the 
DLWC. These electrons hit the FPC close to one of the four WLS fibers 
of the cell. The energy resolution obtained using the RMS and the mean 
value of the distribution {E) is RMS/-E = 8%, while a resolution of 5% is 
obtained from a gaussian fit, which is insensitive to the high energy tail. 



6.3 Tunneling of electrons into the HAC section 

For electrons of less than 100 GeV, more than 99% of their energy is absorbed 
in an infinitely wide block of 26.6 Xq depth. However, electrons hitting the 
WLS fibers start to shower much later than those incident on the lead. 
Since the fibers in the EMC section correspond to only 0.4 Xq, electrons can 
traverse the EMC section through the fibers and deposit their energy in the 
HAC section. For a data sample where electrons are uniformly distributed 
in the transverse plane, about 2% of the events show more than 10% of the 
incident energy in the HAC section. These events are also concentrated near 
the positions of the WLS fibers (see Fig. O). When the FPC is tilted by 



1.7° this fraction of electrons reduces to 0.3% (see Fig. 15). 



7 FPC performance with electrons 
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Figure 13: The left picture shows the distribution of Eppc for 100 GeV 
electrons incident at cell 51. A tail to higher energies can be observed (shaded 
area). In the right picture the positions of electrons as measured by the 
DLWC are shown for the events in the tail. These electrons hit the FPC 
close to one of the four WLS fibers in the cell. 
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Figure 14: The energy deposited in the HAC section of the FPC by 100 GeV 
electrons is plotted for vertical incidence. The plot on the right shows the 
impact postion of electrons depositing more than 10 GeV in the HAC section 
(shaded area on the left plot). 
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Figure 15: The energy deposited in the HAC section of the FPC by 100 GeV 
electrons is plotted for the FPC tilted by 1.7°. 



7.1 Linearity and energy resolution 

For the study of the linearity and energy resolution, electrons incident uni- 
formly on the area of a cell (24 x 24 mm 2 ) have been selected at all available 
beam energies. The total signal has been obtained by summing the signals 
from the cluster of 3 x 3 EMC cells centered on the cell containing the point 
of incidence. 



Figure 16 shows the pulse height distributions for 3x3 EMC cells at different 



beam energies. Gaussian fits have been performed to these distributions. 



Figure |17| shows the gaussian mean versus beam energy and the differences 
from a straight line fit. Two different cell clusters are considered, namely: 
3x3 EMC cells and the whole EMC section. About 91% (99%) of the beam 
energy is observed in the cluster of 3 x 3 EMC cells (EMC section). About 
1% of the electron energy leaks to the HAC section. The deviations from 
linearity are found to be 1% or smaller for both cluster definitions. 

The ratio RMS to mean of the pulse height distribution is plotted versus 



the beam energy in Fig. 18. The data are fitted to the quadratic sum of a 



sampling and a constant term: 
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Figure 16: Pulse height signals from 3x3 EMC cells at different electron 
beam energies. 



RMS _ a 
E ~ VE 



where E is in GeV. The result of the fit is a = (0.410 ± 0.017) GeV 1 / 2 and 
b = 0.062 ± 0.002. 

The signal enhancement around the WLS fiber region distorts the distribu- 



tion and broadens the resolution as discussed in section 6.2. This can be 
demonstrated by selecting only electrons incident on a square of 8 x 8 mm 2 
centered at the cell center. In this way the beam impact point is kept away 
from WLS fibers. Again, the signal obtained from the sum of 3 x 3 EMC 
cells has been considered. The relative standard deviations, cte/E, obtained 
from gaussian fits are plotted versus the beam energy in Fig. [L|. The data 
are fitted to the quadratic sum of a sampling and a constant term: 



E JE 

where E is in GeV, yielding a = (0.34 ± 0.03) GeV 1 / 2 , b = 0.018 ± 0.007. 
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Figure 17: The energy measured for electrons in 3 x 3 cells of the EMC 
section of FPC and in the complete EMC section. The lower plot shows the 
deviations from linearity. 
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By restriction to the 8x8 mm 2 region on the center of the central cell, a 
substantial improvement is obtained in the resolution, particularly in the 
constant term b, which contains the contribution of non-uniformities. 

Lower energy DESY data taken with an FPC prototype are included for 
comparison. There are mainly two contributions to the measured value 
O.SA/y/E: sampling fluctuations and photoelectron fluctuations. Accord- 
ing to MC, the expected sampling contribution to the energy resolution is 
0.29 /y/E. Therefore the contribution from the photoelectron statistics is 
estimated as about O.YI/yE. The LED data taken in those EMC cells, 
which have been used for determining the electron energy resolution, show 
photostatistics fluctuations around 0.15/ y/E, in good agreement with the 
estimate presented above. 



7.2 Position reconstruction 



Position reconstruction for incident electrons is done in the following way: 

• The electromagnetic cell with the maximum energy signal is searched 
for. The event is accepted if around this cell a 3 x 3 EMC cluster can 
be defined inside the FPC. 

• Once the 3x3 cluster is selected, a linear algorithm is used to obtain 
a first approach to the correct position. In the following only the 
position recontruction along the y-direction is shown, but the same 
results are obtained in the x-direction. The formula applied is 



Vba 



J2l=l Si • Vi 
2Ji=l Si 



where i runs over the 9 cells, y% are the y-coordinates of the cells and 
Si are the energy signals. 

The resulting y& ar , when plotted versus the true y-coordinate as given 
by the DLWC, shows a characteristic S-shape as shown in Fig. fl9| . 
After applying a suitable correction to y& ar one obtains y rec that, as 



shown in Fig. IS, gives an unbiased estimate of the true y. The cor- 



rection is of the type 

y rec = Pi+ P 2 ■ tanh(P 3 • y bar ). 
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Figure 18: Electron energy resolution obtained from a gaussian fit to the 
signals measured with a 3 x 3 EMC cell cluster and from the RMS of the signal 
distributions. The beam electrons are distributed either over a 24 x 24 mm 2 
square (cell size) or over a 8 x 8 mm 2 square. The curves show fits to 
the CERN data. For comparison also data from the prototype measured at 
DESY are shown. 
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Figure 19: Position reconstruction for 20 GeV electrons before (upper) and 
after (lower) correcting for the S-shape distortion (see text). 
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Figure 20: Energy dependence of the position resolution for electrons. 



The correction applied at this stage is approximately independent of 
the beam energy, as one could expect from the fact that the shower 



transverse profile is almost energy independent [14]. 



The algorithm for electron position reconstruction explained above has been 
applied to large data samples taken in the CERN and DESY beam tests. 
The position resolution is obtained from a gaussian fit to the Ay distribution, 
where Ay = y rec — yoLWC- The expected 1/VE scaling law is observed, as 



can be seen in Fig. 20 for the y-coordinate. The same behavior is seen for 



x. The result of a fit gives: 

0.82 ± 0.01 



g v = t= © 0.064 ± 0.005 cm 



'E 

where E is measured in units of GeV. 



25 



8 FPC performance with pions 



8.1 Energy response 

The tt~ beam provided by the CERN test beam facility showed a substantial 
muon contamination. These muons have been rejected from the data sample 
by applying the following cuts: 

• Eb 5 < 0.25 mip 

• Eppc > 0.25 • Ebeam 

where E^ is the signal measured in counter 65 (see Fig. ||), 1 mip is the 
most probable signal deposited by a minimum ionizing particle and Eppc 
is the total energy measured in the FPC. 

The 7r~ signal has been measured with different combinations of the FPC 
and FCAL prototype modules. The FCAL prototype modules are 7 A deep. 
For this beam test the modules have been calibrated with electrons and 
muons. Previously, a setup of four FCAL prototype modules had been 
tested with electrons, muons and hadrons, and the calorimeter was found to 
be compensating, e/h = 1.0 for momenta p > 3 GeV to within 3% |J1C[|. 



Figure 21 top shows the signals measured with the FPC and FCAL prototype 
as a function of the ir~ beam energy. For the FPC alone and FPC+FCAL 
prototype the measured signal follows closely a linear rise with the beam 
energy. For FPC+FCAL prototype, the deviations from linearity are <3% 
(see Fig. |2~T| bottom). 

Energy fractions of about 56% (FPC HAC), 75% (FPC: EMC+HAC) and 
93% (FPC + FCAL prototype) are measured at E w = 100 GeV. These 
fractions change slowly with beam energy. 

Signal distributions obtained when combining the signals from FPC and 
FCAL prototype are shown in Fig. ^ at beam energies between 10 and 100 
GeV. 
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Figure 21: T/ie energy measured for ir~ in the HAC section of FPC, in the 
FPC and in the FPC + FCAL prototype modules. The lower plot shows the 
deviations from linearity for the total FPC + FCAL prototype signal. 
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Figure 22: Pulse height signals observed in the FPC + FCAL for ir~ beam 
energies of 10, 20, 40, 60, 80 and 100 GeV. 
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Figure 23: Energy resolution (ge/E) for it - ' as measured by FPC alone and 
FPC + FCAL prototype as a function of the ir~ beam energy. 



8.2 Energy resolution 



Figure 23 shows the energy resolution for tt~ as a function of energy for FPC 
alone and FPC + FCAL prototype. Each point is the result of a gaussian 
fit to the corresponding tt~ distribution shown in Fig. The curve shows 
a fit with the parametrization 



o-e 



E JE 



where E is in GeV, yielding: a = (0.65 ± 0.02) GeV 1 / 2 ; b = 0.06 ± 0.01. 

The resolution is affected by transverse leakage (see Fig. ||). When installed 
in ZEUS the FPC is completely surrounded by FCAL modules. As a result 
there is no transverse leakage except into the beam hole and, according to 
MC, the energy resolution in ZEUS improves by 20%. 
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8.3 Position resolution 



The position reconstruction for ir~ has been studied using three different 
sets of cells corresponding to the FPC-EMC section, the FPC-HAC section 
and the combined EMC + HAC sections. In the following only the position 
recontruction along the y-direction is shown, but the same results are ob- 
tained in the x-direction. A barycenter of the cell coordinates is determined 
as 



Ei w i • Vi 

Ubar 



where i runs over all cells considered, yi is the y-coordinate of the center of 
cell i, and Wi are weights defined as 



Wi = max ^0, In ySi/ Si j - Cij 

Si is the energy signal in cell i and Cj is the threshold in the quantity 
In (5i/Ei Si) corresponding to cell i below which the cell is not included 
in the barycenter calculation. 

The relation between the position y given by the DLWC (yDLwe) and the 
ybar measured in the FPC shows a negligible S-shape distortion. This is due 
to the use of a logarithmic function of the normalized signal to determine 
the weights Wi. Hence no S-shape correction is needed and the reconstructed 
position y rec has b een obtained by just rescaling the barycenter yi, ar . 

The position resolution is obtained from a gaussian fit to the Ay distribution, 
where Ay = y rec — yoLWC- Figure 24 shows the measured resolution, a y , as 



a function of the beam energy. All the cells of the FPC have been used to 
determine the barycenter. The data are fitted using: 



E 



where E is in GeV. The result of the fit yields: a = 22 ± 3 mm- GeV 1 / 2 and 
b = 3.3 db 0.5 mm. When only EMC cells or only HAC cells are used to 
obtain the barycenter, the resolution degrades as shown in Table ||. For the 
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Cells used 


a (mm-GeV 1 / 2 ) 


b (mm) 


EMC + HAC (a) 


22 ±3 


3.3 ±0.5 


HAC (a) 


50 ±4 


2±3 


EMC (E EMC > 0.2 • E beam , RMS) 


74 ±8 


9.6 ±0.9 



Table 3: tt~ position resolution using different sets of cells, with the a from 
a gaussian fit or the RMS, parametrized as the quadratic sum of a sampling 
and a constant term. 



set of EMC cells, only events which deposit more than 20% of the energy in 
the EMC section are considered and the RMS is used instead of the sigma 
from the gaussian fit. 



9 Monte Carlo for the ZEUS environment 



The FPC has been implemented in the Monte Carlo program MOZART 
for the simulation of the ZEUS detector response. The program MOZART 
uses the general purpose packages GEANT 3.13 and GHEISHA modified to 
describe the test beam data measured with the ZEUS uranium calorimeter 
modules. In particular, for the e/h ratio and energy resolution of hadrons, 
discrepancies in the range 20-30% were found between the GEANT 3.13 
+ GHEISHA prediction and the measurements from the FCAL prototype 
[15|. To overcome these discrepancies a shower terminator was introduced in 



the GEANT 3.13 + GHEISHA package, which operates on the evaporation 
energy of excited nuclei and on neutrons with kinetic energies below 50 MeV 



[16|. This shower terminator had free parameters which were chosen such 
that the discrepancies between Monte Carlo prediction and test beam data 
were below 5%. In addition, and because the package GEANT 3.13 was 
very time consuming by tracking shower particles in the calorimeter, a new 
shower terminator was introduced which operates on e~ , e + and 7 with 



kinetic energies below 200 MeV [16]. Using both shower terminators in the 
GEANT 3.13 + GHEISHA package the execution time was reduced by a 
factor of about 15. 

The scheme of shower terminators has also been used to simulate the FPC 
response in GEANT 3.13 ± GHEISHA. Since the FPC and FCAL are differ- 
ent calorimeters, a different set of parameters of the shower terminators has 



been chosen for the shower particles inside the FPC. In Fig. 25 average en 
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Figure 24: The tt~ position resolution, using EMC and HAC section cells, 
as a function of the beam energy. The resolution is calculated as the a 
of a gaussian fit to the difference between reconstructed and true (DLWC) 
position. 
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Figure 25: The mean of the reconstructed energy in FPC, FPC+FCAL, 
FPC-HAC and FPC-EMC for ir~~ is shown as a function of the beam energy. 
The results of the Monte Carlo for the CERN-test geometry are shown in 
comparison with the data. 

ergies deposited in the different sections (FPC-EMC, FPC-HAC, FPC and 
FPC + FCAL prototype) are shown as a function of the beam energy for 
test beam data and MOZART. Good agreement for all sections is observed. 
Energy resolutions (from gaussian fits) are plotted in Fig. || as a function 
of the beam energy for data and MOZART. The FPC energy resolution is 
reasonably well described. When adding the FCAL prototype signal to the 
FPC signal small discrepancies of the order of 5-10% are found. 
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Figure 26: The energy resolution for ir~ in FPC and FPC+FCAL prototype 
is shown as a function of the beam energy. The results of the Monte Carlo 
for the CERN-test geometry are shown in comparison with the data. 
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10 Conclusions 



The Forward Plug Calorimeter (FPC) for the ZEUS detector at the electron- 
proton collider HERA at DESY, has been tested at DESY and CERN with 
beams of electrons, muons and pions in the range 1 to 100 GeV in 1997. 
The results of these tests have been used to determine the calibration of the 
FPC. The energy and position resolutions, oe and a y respectively, measured 
for both electrons and ir~ can be summarized as follows: 

for electrons: a E /E = (0.41 ± 0.02) /y/E © 0.062 ± 0.002 
u y = (0.82 ± 0.01) /yfE © 0.064 ± 0.005 cm 

for pions: o E jE = (0.65 ± 0.02)/ V# © 0.06 ± 0.01% 
G y = (2.2 ± 0.3)/ VI © 0.33 ± 0.05 cm 

where E is measured in units of GeV. 

A monitoring system using a 60 Co source has been used to transport the cal- 
ibration constants from the test beam to the ZEUS detector and to monitor 
the stability of the calibration since its installation in 1998. 
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